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Abstract 



I A dual formulation for the problem of determining absolute performance 

' limitations on overshoot, undershoot, maximum amplitude and fluctuation 

-4— » I 

I minimisation for continuous-time feedback systems is constructed. Determin- 

ing, for example, the minimum possible overshoot attainable by all possible 
stabilising controllers is an optimisation task that cannot be expressed as 
a minimum-norm problem. It is this fact, coupled with the continuous-time 
OO ■ rather than discrete-time formulation, that makes these problems challenging. 

We extend previous results to include more general reference functions, and 
derive new results (in continuous time) on the influence of pole/zero locations 
on achievable time-domain performance. 

1 Introduction 

X 

H ■ In this paper we study the problem of finding limits on the performance of the error- 

response performance, for a specific input, for lumped continuous-time siso systems. 
We are trying to find the best possible tracking performance achievable by a lumped 
(or rational) BiBO-stabilising feedback controller. The theory to be presented is ap- 
plicable to a large class of performance measures of practical significance, including 
overshoot for example. The simultaneous imposition of hard bounds on the output, 
in conjunction with overshoot minimisation, can also be handled, allowing consider- 
ation of the inevitable trade-off between rise-time and overshoot performance. These 
issues have been examined in a discrete-time setting in [5] etc. The investigation of 
time-domain performance limitations is mathematically more challenging in contin- 
uous time than in discrete time. In the continuous-time case it is harder to derive 
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conditions under which the performance hmit for rational controllers is the same as 
for non-rational controllers [rational approximability) . 

The usual approach to these questions, is to express this performance-limit in 
the form of an optimization problem (to be computed via duality methods), over 
an extended ambient signal-space, usually a standard dual Banach space, to ensure 
existence of optimal elements and adequate duahty properties. In general, these 
optimal elements will be non-rational, and hence suffer from implementability issues. 
It is then of importance to find if optimal performance may be approached using 
rational controllers. Since these nonrational elements may be of doubtful practical 
utility, one may well abandon the search for optima, and consider only the class 
of rational controllers in the optimization. This will permit more flexibility in the 
choice of the ambient Banach space since it no longer need be a dual space. 

In earlier works this underlying space was chosen to be (dual of L^), but this 
choice of space does not, of itself, enforce the zero-steady-state-error condition. We 
shall incorporate this condition by using the space Co (of continuous functions of 
time that tend to zero at infinity) as our choice — it does the trick, but it is not a 
dual, so that existence of optimal elements is not assured. This choice of space is 
dictated by the continuity properties of the basic performance measures we study 
here. 

This is the approach followed in this paper. Moreover, to side-step the rational- 
approximation issue, we consider directly the set of error signals due to the ra- 
tional stabilizing controllers, and the optimization over its closure. Under the 
conditions we assume in this paper, passage to this closure has no effect on the 
limit-of-performance. This clarifies the relation between the limit-of-performance 
and rational approximation, but requires an explicit characterization of the closure 
of the feasible set of error-signals. However, it is found that the closure itself has 
pleasant form, well-suited to application of the duality theory. That is, by exploiting 
the structure of this set we have that rational approximability holds by definition. 

Further, we consider a general class of performance functionals and set up a 
general duality framework for the analysis of such problems (in the style of [5] 
for the case of discrete time). These functionals include as special cases the well- 
known classical criteria of overshoot, undershoot and others. This treatment also 
allows a general fixed input (not just a step). The classical criteria just mentioned, 
are continuous with respect to the supremum-norm on the signal space — this fact 
motivating our choice of Cq as the ambient space. 

In the literature to date the optimization is performed over error signals in the 
bigger space L^o, with attendant consideration of rational approximability. In [121 
[T3l [T4] the set of candidate controllers is expanded, moreover, to include those that 
may not be BiBO-stabilising, and the limiting error is unconstrained. Asymptotic 
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stability and zero steady state error are then enforced through selection of suitable 
output weighting. 

Use of the space L^o to formulate the primal has the advantage that the optimal 
solution for the primal problem is guaranteed to exist for the L^o norm minimisation 
problem. For the more general cost functions considered in this paper, however, it 
is clearly not to be expected that an optimal solution will in general exist, neither 
in Loo, nor in Cq. 

We provide a derivation of a dual formulation covering a wide range of problems, 
extending known results. The continuous time norm minimisation problem 
for a fixed input was first considered in [2], but rational approximability was not 
addressed. Miller [7] gave a rational approximation result for response to a step ref- 
erence input. The construction of rational suboptimal controllers for the continuous 
time Loo norm minimisation problem has been considered in [13], [12] and [1]. Yoon 
[9] and [10] extended the class of optimisation criteria, by considering convex combi- 
nations of overshoot and undershoot (with a step-input) and in these latter papers 
is found the first application of conjugate-duality techniques to continuous-time 
control. 

The constraints for the dual formulation presented here can be interpreted as 
arising from an open-loop dynamic system. These structural insights are exploited 
to derive new results on the influence of plant pole/zero locations, or rise-time con- 
straints on achievable performance, giving results of identical form to those obtained 
for discrete time in [5]. The results presented here are applicable for general refer- 
ence inputs ([21 [To] assumes step-input) and to systems with more general exogenous 
inputs, for example a fixed disturbance entering at the plant output. For definiteness 
and simplicity our set-up is the tracking problem shown in Figure [H 

An important design consideration is the extent to which performance objec- 
tives are in conflict. It may not be possible, for example, to find a single controller 
that reduces both the Loo-norm of the error, and the negative-error (overshoot), in 
response to a step, to close to their fundamental limits. Furthermore, the extent 
to which there is conflict depends on the location of the unstable poles and non- 
minimum-phase zeros of the plant. The dual provides such information. Another 
trade-off important in practice is that between rise-time and transient performance 
measures such as overshoot and undershoot. The answer to this question also de- 
pends on the location of the poles and zeros of the plant, and again the interpretation 
of the dual as an open-loop dynamic system provides insight and new results. Also 
in this paper we extend known results on the minimization of the "maximum-peak to 
minimum-peak" value of the error response, that is the difference between the max- 
imum and minimum values of the error response; we term this quantity fluctuation. 
For some applications fluctuation minimization is of more concern than Loo-norm 
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minimization. For some plants there then arises an unavoidable trade-off; either a 
small Loo-norm of the error response, or a low fluctuation, must be sacrificed. To 
what extent they are in conflict depends on the location of the poles and zeros of 
the plant, and again it is the dual that provides answers. 

2 Mathematical Preliminaries 
2.1 basic notation 

We shall write M+ for the real interval [0, +oo). M[s] and M(s) denote respectively, 
the spaces of real polynomials, and real-rational functions, in the complex variable 
s. The set of all proper stable members of M(s) (i.e. those with no pole in the closed 
complex right-half-plane nor at infinity) is denoted by S. The Laplace transform 
of a function /, will be written as /. We let Lp(M+) (1 < p < oo) stand for the 
space of Lebesgue p-integrable functions on ]R_|_. The space of continuous functions 
ip : ]R+ — 7- M for which limt^^ ip{t) = 0, endowed with the supremum norm Hv^Hoo '■= 
supi>o |v5(t)|, shall be denoted by Co(M+). For any a G M, write Co,a(I^+) for the 
subset consisting of those elements ip G Co(M+) for which ip{0) = a. For a subset A 
of a space X, the indicator function la of A is defined on X by 



If / : X M U {+oo}, then dom/ denotes the set {x E X \ f{x) < +oo}. 
2.2 Fenchel duality theorem 

Definition 2.1 Let X be a topological linear space. For any / : X — M U {+oo}, 
the (Young-Fenchel) conjugate f* : X* — )► M U {±00} is defined by: 



It follows that /* is convex and weak* lower-semicontinuous. 

We shall require a Fenchel duality theorem in the following form (see [HI Theorem 
18(a) and Example 11'] for a more general formulation) 

Proposition 2.2 Let / : X — )■ M U {+00} and g : W ^ RU {+00} be convex, 
with X a locally- convex topological vector space. Let A : X ^ M" be bounded linear. 
Assume also that g is finite-valued at some point in A{int dom/). Then, 



and / will be of the form /q + lt for a finite- valued /o and a convex set T. 




for all X* G X* 



f*{x*) := sup ((x, X*) - f{x)) 
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2.2.1 duality for a space of continuous functions 



The application of the Fenchel theorem shall require a characterisation of the dual 
of the Banach space of signals under consideration. In this paper, the ambient space 
of (error) signals will be Co(M+), which has a well-known dual space. Indeed, the 
dual of Co(M+) is iso metrically isomorphic with the space M(]R+) of all regular finite 
signed Borel measures /i on R+, with the variation norm := |/i|(]R4.) (For formal 
definitions of these properties, and a statement of the duality for Co(M+) see, for 
example, [1]). The action oi ^ E M(]R_|_) as a dual functional on Co(M+), is indicated 
by the pairing 



for e G Co(M_|_). (Note also that this expression is well-defined whenever e is bounded 
and (Borel-)measurable, since is a finite measure.) 

3 Problem Formulation 
3.1 primal feasible set 

Consider the set T of all error sequences achievable with a rational stabilising con- 
troller (for the plant P) in the standard one-degree-of-freedom feedback configu- 
ration, see Fig [H Here the plant P is a scalar, linear, proper, finite-dimensional 
system. Mathematically, this condition on P can be expressed by the requirement 
that P G M(s) has a zero at oo of order 6'p > (in the sense that the degree of the 
denominator is the sum of 9p with the degree of the numerator). The reference input 
w is assumed to be an ordinary function w : M+ — > M (no delta-function terms) with 
rational Laplace transform. A typical linear stabilising controller is denoted C. The 
plant output is y, and the error signal is e = w — y 




w + 



o 



e 



u 



y 



c 



p 



Figure 1: A closed-loop control system 



Let n, in iS be a coprime factorization for P. Then d{oo) := lim| 
nonzero and finite, and h has a zero of order 5p at oo, so that 



d{s) is 




lim s ^ h{s) = and lim s ^h{s) is nonzero and finite . 
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Such a factorization for P is readily found; place P = q/r for coprime polynomials 
g, r, with degr = degg + ^p. Then h := q/{- + and d := r/(- + \)'^^^^ both are 

in iS, and by [11, Chapter 2, fact 20], are coprime (in S). 

Given the reference- input then e is a closed-loop error-signal (for some stabi- 
lizing controller C for P) if and only if e = 'wd{v — qh) = w ~ wn{x + qd) for some 
q & S, where x, -0 in iS arise from the coprimeness of h and d, and satisfy 

xh + vd = 1 in S . 

From the assumption on w, its Laplace transform then satisfies 

lim w{s) = , 

Re s— >'+oo 

and w{s) has a zero of order 6*^, > 1 at infinity, implying that lim|s|_^.oo s^'"^^w)(s) = 
and that limi^j^oo ■s^™w(s) is nonzero and finite. Hence, for e as above, e{s) has a 
zero of order at least 6*^ at oo, and from w — e = wh{x + qd) follows that e — w has 
a zero of order at least 9p + 9^ at oo. 

Define our feasible set of possible error signals by the affine set (where the star * 
denotes convolution of functions) 

J-" := {e G Co(M+) \ e = w * d * {v — q * n) ioi some q such that q E S} . 

(We are constraining e to be in C7o(M_|_), as a criterion for zero steady-state tracking 
error). Then, for 9 := {9w,9p), 

T C Xq_ := {e G Co(M+) | e rational, and has a zero 

of order at least 9^ at infinity, and 

e — w has zero at oo of order at least 9p + 9^]}- 

Note that is an affine subspace of Co(M+). By considering the partial-fraction 
expansions of e, it follows that any e G L^(M_|_) U Co(M+) with rational Laplace 
transform may be expressed in the form (for all t > 0) 

e(t) =Re ^Qt'^'e^'* 

i 

where the above sum consists of a finite number of terms, and Cj G C, Re A, < 0, 
ki > 0. 

Introduce the notation for right half-planes 

C+ := {s G C I Res > 0} 
C7:= {s G C I Res > 0} 
c;, :=C;u{oo}. 
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We shall now follow an analogue of the developments of [5]. Let P have poles 
Pi, . . . ,pm and zeros ^i, . . . , z„ in the right-hand plane C+. Also, let the reference- 
input w have zeros Vi,. . . ,vi in C+. It is assumed that none of these lie on the 
imaginary axis {s \ Res = 0}. Further, all these poles/zeros are assumed to be 
mutually distinct, and simple. 
Place 

Me_=leeXe_ e{Pi) = ^ i = 1, 2, . . . , m J> . (1) 





eizi) 


= w{zi) 


2 = 1,2,. 


.,n 




e{Pi) 


= 


2 = 1,2,. 


. , m 






= 


2 = 1,2,. 





Lemma 3.1 With the assumptions as above, T = Mq_. 

Proof: If e G J-", then e G Xo_ as argued earlier. The constraints on e at the pi, Vi 
follow as in the discrete-time case. Conversely, if e G Mg, form q := v/n~e/{hwd) = 
i (-^^ ~ The only possible C+g-poles of q are at the C+-poles/zeros of P, the 
C+-zeros of w, and at infinity. The constraints at each of the C+-points ensure that 
e has no poles there. It remains to check the behavior at oo. Now, d{oo) ^ 0, and 
the prescribed behavior of e and e — w means that = 0{-pg^—^) =0(1) for 
\s\ — 7- OO. Thus q{oo) is finite, and since q has no poles in C+, it is in S. It now 
follows that e G J-". □ 

As in the discrete-time case, whenever Zi and Zj form a conjugate pair, we retain 
only one of these in the list of constraints characterizing Mg. We make a similar 
reduction for the pi and Vi also. This entails no loss of information from Mg (since 
e{z) = e{z) for any z, and any real-valued function e.) 

For each interpolation point zj = Xj + iyj (recall Xj > 0), define 

aj(t) := e'""^^ cosi/jt , aj+i(t) := e"""'^* sin y^t . 

If hj and Cj are also defined similarly with respect to the pj, Vj respectively, Mg 
takes the form 



Mg= (eeXg 



{e,a.i) = {w,ai) i = 1,2, . . . ,n 
{e,hi)= 2 = l,2,...,m }, (2) 
(e,Cj) = 2 = 1,2, . . . ,/ 



where for functions u G L°° and v G L^, {u,v) := v{t)u(t)dt. 
Let A : Co(R+) ^ j^m+n+i defined by 

Ae := ((e, ai), . . . , (e, a„), (e, b^), . . . , (e, b^), (e, Ci), . . . , (e, c,))^ G M™+"+' , (3) 
and let 

6 :=((«;, ai), ... , {w, a„), 0, 0, . . .f e R"^+"+' . (4) 
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Then, Mg has the form {e G Xe_ \ Ae = 1)}. 

Given a performance measure / : — )■ M U {+00}, our question is to evaluate 



(P): 



inf /(e) 



which represents a theoretical limit of performance for "physically realizable" con- 
trollers (in the sense of having rational Laplace transform). This is the central object 
of study in this paper. 



3.2 primal time-domain performance objectives 

We shall principally consider functionals /o on Co(M+) of the form: 



fma ifi 
fposifi 

fosie 

fusifi 



sup |e(t)| = ||e 
t 

sup [(e(t))+] 
t 

sup[(-e(t))+] 
t 



(maximum amplitude) 
(positive error) 
(overshoot) 



sup(e(t)) - inf(e(t)) 



sup [(e(t) — w(t))_ 



(fluctuation) 
(undershoot). 



where for real A, we define A+ := max(A,0) and A_ := min(A,0). It is straight- 
forward to verify that these functionals are all convex, and also continuous (in fact 
Lipschitz, with constant 1) with respect to the || • ||oo^norm on Co(M_|_). 
As in discrete-time, we have fos + fpos = 2//; and that 



fpos 
fos 



< fma < < 2/„ 



Also, for e G Cn 



ffM) = minsup |e(t) - ^| . 



(5) 



(6) 



whose proof follows by trivial modification of the proof of its discrete counterpart 
in [5]. Later, we will consider some simple time-domain constraints, so require 
treatment of functionals of form /o + lt for appropriate choices of sets T representing 
these additional conditions. 

Remark 3.1 Recall that a rational error-signal e{-) satisfies (the zero-steady-state 
condition) limt^oo^it) = iff e E Co(M+) iff e & Lp(M+) (p^ 00). Thus, we could 
have chosen (p ^ 00) as our ambient space — however, the functionals fo above 
are not continuous relative to these -norms, negating the usefulness of this choice. 
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3.3 initial statement of duality for our primal problem 

The analysis of problem (P) will proceed by recasting it in dual form via Propo- 
sition \2.2\ as is given below in Proposition 13.31 As preparation, we require the 
following characterisation of the closure Xg (proved in the Appendix). 

Co_^(M+) {with a = w{0+)) if 0p > 0, = 1 
Proposition 3.2 = { Co(M+) if = 0, = 1 

c7o,o(M+) tf e^>i 

(where Co,q,(M+) C Co(]R+) denotes the set of functions ip for which (^^(O) = a). 

The formulation of the dual problem will also require the following spaces: 

U := span [ai, . . . , a„], V := span [bi, . . . , b„], W := span [ci, . . . , q], 

From the blanket assumptions on poles and zeros, we have f/, V, W contained in 
Co(M+). Also, from the resulting integrability of the a„, b„, c„ , each e* E U ®V(BW 
defines a measure G M(]R_|_) by 

= / e* for Borel sets E . (7) 
Je 

Thus U (BV Q)W may also be considered as a subspace of X* = M(]R+). 
The next basic duality result is the foundation for all subsequent analysis. 

Proposition 3.3 Let T C Co(IR+) be convex, fo '■ Co(M+) M convex and contin- 
uous, with b G A{X0n intT), where A and b are given in ^ and Q). 
If9^ + 9p = l (i.e. Tq_ = CoiR+) via Prop\3lS) then 

inf /o = max [(Proj uil^),uj) - {fo + ixTifJ')] , (8) 
TnXgnA-^b ti&u®vs)W 

and if instead, 9^ + Op > 1 (so that by Prop \3.^ Xg = Co,a(M+) with a = w{0+)) 
then 

inf fo= max [(Proj w) + (/o + iT)*(/i)] , (9) 

TnXenA-^b fj,eU®V®W 

where: Proj u{') denotes the natural projector from U (B V (B W onto U; and, for 
any fi, and any f , 

/#(/i) :=max[aA-/*(/i + A(5)] (10) 
with 5 denoting the Dirac measure at E M_(_ . 



Proof: Now, b G A^Xg D int T) implies, by convexity, that T n XgH A'^b D T fl 
Xg n A~^b. But Xg n A^^ = Xg n A^^b by a simple argument using the finite- 
co dimensionality of A^^b. Thus T fl Xg fl A^^b D T (1 Xg fl A~^b and hence by 
continuity of fo, 

inf fo = Jnf fo . 
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We now apply Proposition 12.21 to the latter minimization in Xq_. Place f '■= fo + lt 
and g := If = Cq, and taking X = Cq, then b G yl(int dom/) and after 

some simple manipulations of the resulting dual, we obtain the form ([8]). 

If instead, JQ = Co,a, note that A-'^b n X^ = A~% where Ae := (Ae, e{0))^ (for 
all e G Co(M+)) and b := (6, a)^. We also have b G A(intT), so we may apply the 
duality to infrnCo.cnA-ib /o = ^^^rnA-H /o, with X = Co again, with the new A, b, 
which eventually yields the dual in (jH]). □ 

4 Dual Formulation 

The dual characterisation will be completed by evaluation of the conjugate function- 
al appearing on the right-hand-side of equations ([8]) and ([9]) above. We now study 
the forms of the conjugates for a general class of cost-functions which includes those 
of interest in this paper. 

We assume the objective functional / : Co(M+) — ^ M U {+cx3} take the form 

/(e) = supF,(e(t)) (11) 

where F satisfies the following assumptions. 

1. For each t > 0, Ft : R ^ [0, +oo] 

2. For all t > 0, all L > 0, the sublevel-set [Ft < L] is a nonempty closed (possibly 
unbounded) interval in M, with endpoints e^{L) G M U {— oo}, e^{L) G M U 
{+00} respectively, with, further, 

infe^(L) > —00, supe7(L) < +00. 
*>o t>0 

3. For each L > 0, the functions t 1— )■ ef{L) are piecewise continuous (in appro- 
priate sense for extended-real-valued functions) with at worst a countable set 
of jump- discontinuities, where at all such jumps t, have 

min{e^:!:^(L), e'^^L)} > max{ef^(L), e^jL)} . 

4. For each L > 0, there is to^L) > such that for all t > to, have that et{L) < 
< e+(L). 

Before we proceed further, a review of some relevant measure theory is appropri- 
ate. (Our reference shall be [1].) Given /i G M(M_|_), we can find a Hahn decomposi- 
tion of M+ into disjoint Borel sets P, N such that n{E fl P) > and /i(E fl X) < 
for each Borel set E. From this are obtained regular finite Borel measures ii± given 
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by n+{E) := n{E H P) and n^{E) := fi{E D N), with fi+ > and < 0. We then 
have fi = fi+ + fi- and the variation of fi (denoted \fi\) is defined by := /U+ — 
Note that our sign convention has been chosen to conform with that used in the 
discrete-time analysis of |5j, but differs from the usual choice in measure theory, 
where and /i_ are both non-negative, whereas here we have /i_ < 0. Recall that 
we assume M(]R+) to be normed by ||/i|| := |/i|(M+). 

Note that if measure E U Q)V Q)W and function e* are related via (JTj), then by 
standard arguments follows that 




for Borel sets E , 



and that /i > if and only if e*{t) > for all t, with a similar relation for the reverse 
inequality yU < 0. 

We also recall the following standard measure-theoretic convention: Since we 
need to integrate real functions taking +oo as a possible value, the definition of the 
integral (with respect to a measure) incorporates the convention 0-cxd := whenever 
one of the factors is the value of a (unsigned) measure. 

4.1 The Conjugate for the Assumed Form of the Primal 
Objective 

The following result provides a continuous-time analogue of part of [5l Theorem 9] . 
(The proof is deferred to the Appendix.) 

Theorem 4.1 For f : Co(M+) — R U {+00} defined as in ( ITTi) . where F satisfies 
the assumptions 1 — ^ above. Then, 

/ et{L)df,^{t)+ e;{L)dfi.{t)-L) (12) 

for all G M(]R_|_) such that = for each t for which there is some L for 

which at least one of e^{L) or e~'~(L) is discontinuous at t. 

Remark 4.1 If Ft = F (independent oft), assumption 3. is inactive, and as- 
sumption 4- can only be valid for L satisfying L > F{0) (which is equivalent to 
[/ < L] 7^ 0^ and from the proof (see Appendix) it follows that the supremum in 
( fl^) is then to be taken over L > F{0). 

4.2 Duals for some Time-Domain Minimization Problems 

We may now derive the duals of the continuous-time versions of MA, OS, POS, FL 
in the manner of [Sj. 
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Because of Proposition I3.2[ this separates into cases where Xe = Co(M+) and 
= Co,a(IR+)- We begin with the former case. 

4.2.1 Maximum Amplitude 

Clearly e^{L) = L and e~{L) = —L. By Theorem 14.11 and Remark 14. H 
dom/;, = G M(M+) I M < 1} and fUf^) = 0. 



Thus by Proposition 13. 3[ the dual of the problem of maximum-amplitude mini- 
mization, denoted MA DUAL, can be written as 

max (Proj (/(/i), w). (13) 

ti&U(BV®W 

IImII < 1 

4.2.2 Positive Error 

For this case we have e~{L) = —oo and e+(L) = L. By Theorem l4.1l and Remark |4.H 
dom/;„, = {/i e M(R+) I /i > and /i(M+) < 1} and = 0- The dual of POS 

is 

max (Proj (14) 

fi&U<S)V(BW 

M > 0, IImII < 1 



4.2.3 Overshoot 

For overshoot minimization e~^{L) = +oo and e~{L) = — L. By Theorem 14.11 (and 
the Remark), dom/;, = {/i G M(M+) | /i < and fi{R+) > -1} and f*^{fi) = 0. 
The dual of OS is 

max (Proj (/(/i), w). (15) 
At e c/® ve w 

< 0, < 1 



4.2.4 Fluctuation 

Also, as in the discrete-time case, we may use ([6]) to deduce that f^i is the indicator 
function of 

{/i G M(M+) I /i+(R+) < ^, /x_(M+) > , 
(see proof of O Theorem 11]). Hence, by Proposition 13. 3[ the dual of FL is 



max (Proj [/(/i), w) . 

iieu®vs)W 

A'+(R+) < I 
At_(IR+) > - i 
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We observe in passing that for the fluctuation-minimization problem, that the 
minimum is never achieved (even over Mg), except in the trivial case where P has 
no poles in C+, so G M := Mg is optimal. Indeed, suppose an optimal e G M 
is attained, and let fi denote the optimal dual element. These must satisfy the 
alignment condition ffi{e) + fji{fi) = {e,fi), which yields 

(e, /i) = ffi (e) = ^ (sup e - inf e) 

since fjiifi) = 0. If FLopt = this immediately yields the contradiction e = (since 
^ M). If FLopt is positive, note that the optimal /i is then nonzero. Now as 
/ dfi+ < I and / dfi- > — |, and sup e > > inf e, we obtain 

> sup e J d^+ + inf e j dfX- , 

so that 

< y ^ ^) '^/^+ ^ J ~ '^''^^ — ^ ■ 

Since both terms in the above are nonpositive, we conclude that 
y (e — sup e) dfx^ ~ J ~ ~ ^ ' 

Writing this in terms of the function e* associated to fi via (^^), we have for almost 
all t (w.r.t. Lebesgue measure), and hence by continuity, for all t, that 

(e(t) - supe)e*(t) = = (e(t) - inf e)e*(t) . 

Thus 

j sup e if e* (t) > 
^^^^ ~ [ inf e if e*(t) < 

By analyticity of e* 7^ (being a finite sum of sinusoids), each of its zeros is isolated, 
so by continuity of e follows that inf e = sup e, implying again the contradiction 
e = 0. 



4.2.5 Undershoot 

We assume w{t) > for all t. (Note: To derive an expression for f*^, we do not 
yet need to assume a rational Laplace transform for w). Place Ft{^) := (^ — w(t))^. 
Then Ft{0) = for all t > 0, and G {Ft < L} = {-00, w{t) + L] for each L > 0. 
Theorem 14.11 gives, for all fi G M(]R+), that 

f:M= I w^f//i+ + (-oo)-/i_(M+)+supL(/i+(M+)-l) 

which yields the following 
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Corollary 4.2 Let w>0. Then, 

dom/:,C{;,GM( 

and for each /x > 0, 



J w dfi if fi > and ||/i|| < 1 
+00 otherwise 



Note that j wdfi may take an infinite value. 
Consequently, the dual for US has the form 



max {PToiv(Bw{fJ'),w). 

< 0, < 1 



(16) 



Recall that each fi E U Q) V (B W corresponds to a function e* by Equation ([7]). 
This has the consequences: \\fj,\\ = |/i|(]R+) = / |e*| = ||e*||i; /i > (as a measure) if 
and only if e*(t) > for all t; and /i < if and only if e*{t) < for all t. We may 
therefore restate these duals in a form identical to those in discrete time: 



(MA DUAL) 
(PCS DUAL) 
(OS DUAL) 
(US DUAL) 
(FL DUAL) 



4.2.6 case of Xa 



Co,a 



max (Proj [/(e*), w) 

e* GU ®V ®W 
\\e*\\i < 1 

max (Proj u{e*),w) 

e* > 0, lle'lli < 1 

max (Proj u{e*),w) 

e* < 0, lle'lli < 1 

max (Projyevy(e*),w) 

e* eU @V (BW 
e* < 0, |le*||i < 1 

max (Proj (7(e*), w) . 

e* eU (BV ®W 

/r+ - 2 



When Xq = Co,q, we apply the duality formula @, which requires calculation of 
/*, given by (fTOll . The forms of all the duals, except for FL, will be unchanged, 
because of the next result, proved in the Appendix. 
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Lemma 4.3 /// stands for any of the junctionals treated (except for ffi), then for 
fi e U ®V ®W, have 

/#(/i) = -r(/i), 

whereas for FL, we have 

ffM = -rjM + «(i/2 - • □ 

This, for FL, yields the dual 

(FL DUAL) max (Proj (/(e*), u;) + a(l/2 - / e^) . 

e* eu ®V ®W Jn+ 

/r+ - 2 
/m^. — ~2 

4.3 Overshoot/Undershoot Minimization 

We follow the procedure of 0. The first lemma is essentially the continuous-time 
counterpart of [5l Proposition 15]. 

Lemma 4.4 Let yi, . . . ,yN 7^ with yi 7^ ±.yj whenever i ^ j ■ If ai, . . . ,aN and 

Pi, ... , are real scalars such that for some C, p > we have 

N 

^^(oj COS yit + (3i sin yit) > —Ce^^* for all t large 

i=l 

then ai = = (3i for i = 1, . . . , N . 

Proof: Place a{t) := Xlili cos + /3jsin?/jt). For r > 0, set 



S{r) 



POO N 

/ a(t)e-''dt = y^ 

Jo t:f ^ 



UiV + (5iyi 



Since Vi ^ for all i, there exists A > such that |5'(r)| < A for all r > 0. Thus for 
some to and any N eM. and r > 0, 

r-JV 



0< / a+{t)e-''dt< / a+{t)e-'''dt 
Jo Jo 



00 



S{r) - I a_{t)e-''dt 

"to 







< A + / \a^{t)\dt+ / \a-{t)\dt 

Jo J to 

/•to roo 

< A+ / \a^{t)\dt + C / e-P^dt 

Jo Jo 



to 



<A + / \a.{t)\dt + C/p:= \' 
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For each fixed A^, we may let r — )■ in the above, to obtain by the Dominated 
Convergence Theorem that < aj^(t)dt < A'. Since |a-(t)| < Ce^^^ for all large 
t, it follows that \a{t)\dt is bounded above for all N, so that a(-) is in L^. This 
implies that the Laplace transform d{s) is defined for all s such that Res > 0, and 
is bounded in this region. However, for all s with positive real part. 



ais 



and so if (a^, f3k) (0, 0), then from the assumed distinctness of the i/i, follows that 
d must have a pole at s = ±iyk, contradicting the boundedness of the transform at 
points of the imaginary axis. Thus necessarily {ak, (3k) = (0,0). □ 

Lemma [4.41 can now yield an analogue of O Proposition 17] in continuous time. 
Given the complex numbers and fi,.., let 7 denote the smallest value 

among those poles/zeros that lie on the positive real axis; if there are none, set 7 to 
00. 

Theorem 4.5 Consider the 'equivalence' class A, of all plants which: 

1. have the same 7 value; and 

2. have the same poles and zeros in the region {s | Res > 7}. 

Then, for any fixed minimum-phase reference input (i.e. having no zeros in the 
closed right half-plane), all plants in A have the same value for OSopt- If the 
reference-input is also non-negative, then all plants in A have the same value for 
USopt ■ 

Proof: We show that the poles/zeros in {s | Res < 7} do not affect the value of 
the dual. 

For simplicity, we assume that all the poles/zeros in {s | Res < 7} have equal 
real part. The general case is handled by iteration of the argument given below. Let 
e* eU ®V ®W with e* > (or e* < 0). Regrouping the modes of e* according to 
the size of the real parts of the associated poles and zeros gives 

e*(t) = e-^^* 5^(af ^ sin(yWt)+/3« cos{yfh))+e-^-' Y.^af^ sin(yf t)+/3f ^ cos{yfh))+. 

i=l i=l 

where < pi < p2 < and in the first sum, which corresponds to the poles/zeros 
in {s I Res < 7}, we have yf'^ 7^ 0, and yf"* 7^ '^vf^ whenever i 7^ j. Hence 

J](af ) smiyi'h) + /3f ^ cos(yf ^t)) = e*(t)/e-^^* + e^^-''^)* ^ . . . 

i=l i=l 

> + 0(e'^P'-P^^') . 
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By Lemma |4.4[ al = = I3\ for i = l,..,Ni. Thus the modes associated with 
poles or zeros with real part smaller than 7 do not contribute to the dual, as claimed. 
□ 



4.4 Analytical Results for a First— Order Plant 

We first derive analytical expressions for performance limitations in terms of pole/zero 
locations for first-order plants. Thanks to Theorem 14. 5 [ in the case of overshoot and 
undershoot limitations on performance, these analytical results can be extended to 
include cases where the plant has an arbitrary number of oscillatory poles or zeros 
(i.e. those off the real axis), and the reference input has an arbitrary number of 
oscillatory zeros. 

Proposition 4.6 Suppose the plant, P = jE^, has one real positive zero, zi^ and 
one real positive pole, pi, where z\ > pi. Let w{s) = 1/s. Define h := pi/{zi — pi). 
Then 

(i) OSopt = h 
(iz) MAopt - ^ 



{Hi) FLopt 



1 - 



2h^ 

(tv) POSopt = 1. 

Remark 4.2 It can be verified using elementary calculus that, for h G (0, 00) 
max(l, h) < -— ^ < ^ r-i < 



1 _ 2-1//^ - h'' - 1 - 2-^/^ 

From (a) and (Hi) it follows that 2MAopt = FLopt if o.nd only if h = 1, that is 
Z\ = 2pi. If h ^ 1, then FLopt is strictly less than 2MAopt; in this example, when 
h ^ 1, a minimal fluctuation response is obtainable only by allowing the infinity- 
norm to be larger than optimal. 

Proof: We have A = span [e~^^*] , B = span [e^^^*] , C is empty and w = 1 (the 
unit step function). 

(i) The dual optimal vector for (OS DUAL), e* = ({t) — r]{t) where C ^ ^ ^-^id 
rj E B, will satisfy ({0) = f]{0) , so ({t) = ae~^^* for some real number a and 
rj{t) = ae~^^*. Hence a j^{e~'^^^ — e'~^^^)dt < 1, from which a < piZi/{zi — pi). 
Then OSopt = niax^g^(^, r) = max^ C{t)ldt = maxa [a/zi] = h. 
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(ii) (MA DUAL) is 



max [a/^i] 

roo 

subject to / |ae~^i* - /Se^^^*! dt<l. 
Jo 

It is clear that, at optimality, a > (3 > 0, and there will exist a positive number to 
such that C(^o) — v{to)- It is obvious also that at optimality the constraint inequality 
will be satisfied as an equahty. The dual becomes 

max 

/•to poo 

s.t. ae-^i* = ^e-^'i* and / {ae'''^ - pe'P'^) dt + {^e-P'^ - ae'''^) dt ^ 1. 

Jo Jto 

This is a simple finite-dimensional constrained optimization problem, which can 
be solved using the method of Lagrange multipliers. After some elementary algebra 
one obtains to = ^^J^ = = 2'^ and a = zi/[l - 2(1/2)^i/pi], from 

which the result follows. 

(iii) At optimality the inequalities in the two constraint equations for (FL DUAL) 
will be satisfied as equahties. Then (FL DUAL) becomes 

max 

poo poo 

s.t. / \ae-''* - pe-P'^\dt^l and / {ae'^'^ - pe-P'^) dt ^ 0. 
Jo Jo 

The second constraint gives immediately that api = /3zi. Let to be the positive 
number with the property that C(^o) = ^(^o)- After performing the integration in 
the first constraint, and writing it as an equality, the dual becomes 

max [a/zi] 

a, to 

S.t. 1 = — (e-P''° - e-'''°) . 

Zi 

After some algebra the optimizing to and a are found to be to — 
a = Zi{l + l/h)^{l + h)/2. 

(iv) In (POS DUAL), if r > 0, an optimizing t] will be identically zero. Then the 
optimizing ( will be positive and satisfy || ( ||i= 1. For w =1 the optimal dual cost 
is cut = 1. □ 
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4.5 A time-domain constraint effect 



We consider some relationships between overshoot (and undershoot) with some sim- 
ple finite-time-horizon constraints. (We have in mind the effect of rise-time con- 
straints on optimal overshoot performance.) Time-domain signal bounds will be 
represented by the set 

T := {e G Co{R+) \ 0-(t) < e{t) < 0+(t) Vt > 0} (17) 

for suitable bounding functions : IR+ — t- M U {— oo,+oo}. For simplicity, we 
restrict to the finite-horizon case, where for some positive t > 0, have (j)±{t) = ±oo 
for all t > i. Also, assume that (j)± are finite- valued and continuous on [0,f|, with 
4>-{t) < <f)+{t) for all t, and 0+(t) > 0. In this case, T has the form 

T = {e G Co(M+) I 0_(t) < e{t) < 0+(t) Vt < t} . (18) 

By selecting a specific e G Xg such that Ae = b, and selecting (p±(t) such that 
4>-{t) < e{t) < (j)+{t) for all t <i, we may ensure that 

be AiX^nintT) (19) 

which suffices for duality to hold, if the cost-function is of the form /o + lt (where 
/o is any finite- valued functional, such as one of those listed in Section [3?2l) . 

Suppose also that the reference signal w is minimum-phase, and that all the 
unstable poles and zeros of the plant are oscillatory (that is, they all have non- 
vanishing imaginary part). As observed in Theorem 14.51 when T is absent, have 
OSopt = 0. Now introduce the further constraint represented by the set T, and 
consider 

OSL := inf /os(e)=_inf (fos + i^T)- 
eeTnXepA-H XenA-'^b 

(where the latter equality follows from flTU]) ). With reference to Proposition 13. 2^ 
make the additional assumption in the case when Xg = Cq^q,, that a > 0, in which 
case we then have 

</>+(0) > a> max{</)_(0),0}. 
(remember that in this case a = w{0+) = e(0)) 

Proposition 4.7 With the above assumptions, OS^p^. = 0. 

Thus, in particular, for step-input, the imposition of rise-time constraints does not 
degrade the optimal overshoot performance. 

Proof: Now, /(e) is of the form sup,>o Ft{e{t)), where F^i^ = i^°^(0+^ (<^_(t),<A^w) (0 
from which follows (since 0+(t) > 0) that ef{L) = (p+it) and e^(L) = max{0_(t), —L}. 
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and that Assumptions 1 — 4 for Theorem 14.11 are satisfied (note, for assumption 2, 
needed 4>+(t) > to get [Ft < L] 7^ 0), and hence for any fi e U ®V ®W, have 

= / + ( + CX)) ■/i+((t, +00)) + 



+ sup 

L>0 



max{</)_, —L}dfi^ — L{fi^{(t, +00)) + 1) 



> / 0+ + / max{0„, 0}(i/i_ + (+00) ■ +00)) 
Jo Jo 

where the first two terms of the latter are clearly finite. Thus, whenever /i G dom /* 
it must follow that +00)) = 0, and so dom/* C {^u | fi\(t,+oo) < 0}. 

The dual then has the form (when Xe = Co) 

max [{PToju{e*),w)-ne*)], 

e- £U ®V (BW 
e* € dom /* 
e*|(t,+oo) < 



but each dual e* now satisfies e* = by Lemma 14. 4[ by the same argument as used 
in Theorem 14.51 Hence, since /*(0) = sup2,>o —L = 0, the dual has value 0. 
If, instead, Xg_ = Co,a, the dual takes the form 

(D) = max [{PToju{e*),w) + f*{e*)], 

e* & U eV (BW 
e* G dom(-/#) 

where 

/*(/i) = max(aA - f*{ii + \5)) . 



We want to show that dom(— /*) C {ji \ /i|(f^+oo) < 0}. Now, we may apply 
Theorem O for /*(/i) for any e U ®V ®W ®M.5 (since for any L > 0, t f-> ef{L) 
has no discontinuities at t = — only at t = t > 0, where |/i|({t}) = 0). Thus, 



-/#(/i) = min(-aA + r(/i + A5)) 

> min ^ — aA + j (j)+d{jji + \5)+ + j max{0_, 0}(i(/i + A5)_ + 
(+00) ■ (/i + A5)+((t,+oo)) 



[+CX)) ■ yU+((t, +00)) + / + / max{0„, 0}(iyU_ + 

Jo Jo 

+ min(-aA + A+0+(O) + A_max{0_(O),O}) (by Lemma E^]) 



agk 



(+00) • yU+((t, +00)) + / + / max{0„, 0}(i/i_ 

Jo Jo 
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where we used the relation </)+(0) > a> max{</)_(0), 0} to ensure the latter mini- 
mization has value zero. 

Therefore, if —^^{fi) < +00, then +00)) = so that /i|(t,+oo) < 0, which 

implies, since fi E U Q) V Q) W, that fi = 0, whence (D) = /*(0) = maxAgM(aA — 
n\6)). But, n\6)) = sup^>o(A+eo+(i:) - L) = A+0+(O) since e+(L) = 0+(O) for 
all L, yielding (D)=maxAeR(aA — A+0+(O)) = as « < 0+(O) (and a > 0). Thus, 
OSopt = in this case also. □ 

5 Conclusions 

Using a dual formulation, new results on fundamental time-domain performance 
limitations for continuous-time systems have been presented. For the problem of 
designing a feedback system to optimally track a specific input, or reject a specified 
disturbance, there are many time-domain performance measures of the output signal 
that can be used. In addition to overshoot, undershoot and the infinity norm of the 
error signal, a performance measure of practical significance, termed fluctuation, has 
been investigated for the first time in a continuous-time setting. 



6 Appendix A: Proof of Theorem 14.1 



Firstly, we consider the case where, for each L, the ef{L) are bounded in t > 0. 
Now, = sup2,>Q supggc-ij j-(g-)<j;^((/i, e) — L) , noting that if F is independent of 

t, the supremum is restricted to L satisfying L > -F(O), as indicated in Remark 14. 1[ 
We fix a value of L and evaluate the inner supremum, and will show that 

p+oo r+oo 

sup (/i,e)= / e+(L)rf/i+(t) + / e~{L)dix_{t) 
e e Co "'0 Jo 

/(e) < L 

from which (|T2|) then follows. 

For fi G M(]R_|_), let {P,N) be a Hahn decomposition [1] of M+ relative to the 
measure /i. Now, among all e(-) : IR+ M satisfying e(t) G [e^ (L) , ef (L)] for all t, 
{n, e) is maximized at e, where 

e(t) := 



e+(L) teP 
ei{L) teN 



This maximal value (/i, e) at majorizes the required quantity. Note that e is bounded, 
so J \e\ d\fi\ < oo, but generally not continuous, nor does it decay as t — )■ -|-oo. It 
remains to show that this upper bound is approached for e G Cq. 
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Let e > 0. By regularity [T] of the measure |e^ — e \dfi, there are Kp C P C Up 
and Kn C C Un {U open, K compact) for which 



Suitable shrinkage of Up, Un can ensure that Up fl Kn = ^ = Un ^ Kp (since 
we can take P and to be disjoint). Let {{pp,ipN} be a partition-of-unity (see 
[3]) relative to the open covering {Up,Un} for R+. We could now propose that 
e{t) := {L)ipp{t) + e^{L)ipNit), so that e{t) G [e7,e^] for all t. However, e again 
may not be continuous. To remedy this, we seek suitable continuous approximations 
for e"*" and e~. 

Let t > be any point of (jump-) discontinuity for (e"*" or or both). For 
ease of presentation, we assume strict inequality in assumption 3. 

For definiteness, assume e^_ > ef^. Form a line-segment, of large negative slope 
—A, from (t,e^^) G to the point of first contact with the graph of e~^\[Q^t). Using 
this, we define a function p+, whose epigrapl^ is the union of the epigraph of e"*" with 
the set of points above the indicated line-segment. This yields a function continuous 
at t, satisfying p"*" < e"*" on R_|_, with p"*" = e"*" except on a small interval [i — a{X),t\ 
where a(A) — > as A — t- +oo. Further, since > e^^, on taking A large enough, 
we ensure, by the continuity of e7 near (but not at) i, that also p'^it) > e^, for all 
t. An analogous construction applies (using line-segments of large positive slope) if 
instead < ef^. 

Similarly, we obtain continuous at t, for which e~ < p^ < on M_|_ and 
p^ = e~ off a small interval at i. It now follows that 



Repeating this process at each such t (these are countable in number) yields con- 
tinuous functions on M_|_, again denoted by p^, such that e~ < p^ < e"*", with 



[Remark: if instead, have equality in Assumption 3. then the above line-segment 
construction will not do — in this case one may easily adapt this, by using nonlinear 
segments that are local graphs of continuous functions.] 

Now, we may define continuous functions e, by e := p'^{t)ipp{t) + p~ {t)ipNif)- 

^the epigraph of a function f : X [— oo, +oo] is defined to be the set {(a;, A) G X x M | f{x) < 
A} of points 'above' the graph of /. 





p — e I d\p\ = constant ■ [jump at t for e ] ■ |p|({0}) = . 



/|p±-e±M|p|<e/2. 



22 



Since the < ipp^N < 1 and ipp + = 1 we get that e < e < . Also, 
|(^,e) - (/x,e)| < y |e-e|(i|/x| 

< / le"*" — e~| + / Ip"*" — e"*"! + / \p~ — e~\ d\iJi\ 

J{Up\Kp)LI{Un\Kn) Jkp Jkn 

< 2e 

At this stage, the proof is not quite complete, since e may not tend to as t — )■ oo. 
For the final step we form a sequence of truncations e„ such that: e„(t) ~ for all 
t > n+1; they are continuous; and satisfy {/i, e„) — )■ (/x, e) as n — >■ oo. We will ensure 
continuity for these truncations by again using steep interpolating line-segments. 

If e{n + 1) = 0, we merely take e'„(i) — ior t > n + 1. If instead, e{n + 1) 
is nonzero, we join the point {n + 1,0) e to the graph of e|[o,n+i) by a steep- 
but-nonvertical line segment, to yield a function that vanishes on [n + 1, +oo) 
and agrees with e on [0,n], with |e„| < |e| on IR+. Also, since e [F^ < L] for all 
t > to, then on taking n > and noting that e{t) G [Ft < L] for all t, we have 
e{t) e conv{0,e(i)} C [Ft < L] for all t. Thus, e„ e Co(K+), /(e„) < L, and 

3„ — e| = / |e| + / |e„ — e| 



< 3 / |e| since |en| < |e| 

7 n 

^0 as n — )• OO , 

and so it follows that 

(a*, e) < (a*, e) + 2e < (/x, e„) + 3e < sup (/x, e) + 3e . 

e e Co 
/(e) < L 

As e > is arbitrary, this establishes the required equality, and hence the formula 
for for the case of bounded e^. 

For the unbounded case, for each R > 0, let Ff- :— Ft + L[-r^+r], where the latter 
denotes the indicator-function for an interval. Then the sublevel-sets of F^ have 
the form [F^ < L] = [ef-(L), ef+(L)] where 

ef+(L) := min{i?, e+(L)} , and ef-(L) := max{-i?, et{L)} . 

Define /i?(e) := supj>Q F^{e{t)). Then, for any allowed fi (since the e^^ are bounded 
and fa satisfies Assumptions 1 — 4) we may apply the previous argument to Jr, to 
yield 

(r+oo r+oo \ 

/ ef+(L)d/.+ (i) + / ef-(L)d/._(i)-L 
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for any allowed ^. By an easy check, for any e G Co and any R > sup |e|, /(e) 
/^(e) = infij>o fnie) = limR^^ fn{e) (since fn I as i? f) so that 

(p+oo r+oo 
/ e^~^{L)dfi+ + / e^~{L)dfi- — L 
Jo Jo 

(r+oo r+oo \ 

/ e^{L)dfi^ + / e~{L)d^_ — L ) 
Jo Jo J 



the latter equality (in MUj+oo}) following from the monotone convergence e^'^ ^ 
and I e~ . 



7 Appendix B: Proof of Proposition 13.2 



The proof will require a couple of preparatory lemmas. We remind the reader of the 
following Notation: 

R+ := [0,oo) 

Co(]R+) := {y? : IR+ — 7- M I (y9 continuous, f{c>o) = 0} 

Co,a(M+) := {if G Co{R+) I ^(0) = a} 

9 := {9p, 9^), where 9p > 0, 9^ > 1; and 9 ■=9p + 9^ 

It will be seen that the closure of Xg is always one of Cq or Co,a (for suitable a), 
depending on the value of 9^,9 p. 

We start with a simple result for Laplace integrals. 

Lemma 7.1 Le^ : M+ ^ M such that Lp'-''\t) = 0{e''^) (for some a > 0) and 
^9^*^(0+) exists and finite, for i = 0,1, . . . k. Then for all s E C for which Re s > a, 

^(s) = ^{s)s'' - ip{0+)s''-^ - (^'(0+)s^-2 - ... - if^''~^\0+)s - ip^'"^\0+) . 

Proof: Integrating by parts yields (p^^\s) = s(p^^~^'>{s) — ip^''^^\0+) for such s. 
Repeat for yj'^'^"^) etc. □ 

We assume that the function w has rational Laplace transform. This ensures that 
m;W(s) -> as Re s -> +oo, and w^'\0+) = w^'\0), for alH = 0, 1, 2 . . .. Moreover, 
Lemma mi applies, for every k. 

Now, define 

Xq := {e G Co(]R+) | e rational, e(oo) = (order at least 9^) 

(e — w){oo) = (order at least 9 = 9^, + 9p)} . 
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Also, form the subspace 
Ye_:= |e:t^^Re^.g,Qt^'e-^'* 
Lemma 7.2 ^ X^. 



Cj G C, fcj > 6'iu — 1, Re Aj > 0, / finite, 
eO)(0) =wO)(0+), j =0,1,...,^ -2 



Proof: If e G l£, then clearly e has a zero of order at infinity. Now, [(e — 
— 7- when Re s — )■ oo. By Lemma [711 applied to := e — w, we obtain. 



w 



when s — )■ oo, 

^ [(e - W7)('~^)ns) = s'"'(e(s) - w{s)) - /"'(e - w)(0+) - ... - (e - m;)('-')(0+) 

= /-^(e(s) 

on using the constraints at t = 0+. Hence e — w has a ^-order zero at oo. □ 

Lemma 7.3 If e e Xg, then: e^^\0) = for j = 0, . . . , 0^- 2; and also e^^\0) = 
w(^)(0+) for J = 0,..., 9 -2. 

Proof: Let k := 9 — 2. Now [(e — — )• at infinity. From Lemma mi 
sfc+i(e(s)-u;(s))-(e(0)-^(0+))s'=-. . .-(e('=)(0)-w('=)(0+)) = [(e-wY'+^^ris) 
0. By assumption on e, have s^^^(e — w){s) — at infinity, and this forces the 
equalities e(0) = w {0 +),... ,e^''\0) = w^''\0+). For j = 0, ... ,9^-2, argue as 
above using e instead of e — ly, to get e*^''^(0) = 0. □ 



Remark 7.1 If9w > 2, then Lemma \77^ requires ^^■'•'(0+) = for j = 0, ... ,9^ — 2, 



otherwise the statement of the lemma entails a contradiction, and Xg will be empty, 
hence of no interest. If 9^ = I, the conditions e^^\0) = (j < 9^ — 2) are absent, 
so w(0+) is unrestricted in this case. 

For any ^, form the set We_ ^ by removing the constraints at t = in the 
latter, that is. 



We_ := I e{t) = Re ^ Qt'^'e"^^* 



Cj G C, /cj > 6*^0 — 1, Re Aj > 0, / finite > . 



Observe that e{t) = 0{t'^) at t = (where m := 9^-1), it follows that e(^)(0) = 
for j = 0, 1, . . . , m — 1 when 9^ > 1. 

We turn now to the proof of the Proposition. Because of the preceding lemmas, 
it suffices to prove density for the subsets Yq of Xg. 
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7.0.1 Case of 6^ > 1 



If 9uj > 1, note that ^^(O) = = {p{+oo) whenever ip G Wq. Also, We forms a 
algebra of continuous functions on the one-point compactification (see [3]) [0, +oo] 
of since {p{+oo) = limj_j.oo exists for each such (p. We seek to apply the 
Stone- Weierstrass Theorem (see, for instance, [3], or any functional analysis text) 
to establish the density of We in a suitable subspace of C([0, +oo]). However, there 
is a difficulty, in that We does not 'separate all points' of [0, +oo] in the sense of the 
Theorem — indeed, by definition of We, the points and +oo cannot be so separated 
(if (p G Wg, it attains equal values at and at +oo). This is the only pair of points 
that cannot be separated in this sense. We resolve this issue by insisting that and 
+00 are 'the same point'. 

To do this, let X denote the the quotient space, formed by identifying the points 
and +00 in [0, +oo], and endowing this with the quotient topology |3j. Then X 
is also compact, being the image of [0, +oo] under the quotient map: [0, +oo] — X, 
which is continuous by definition of the quotient topology. It is easily verified that X 
is also Hausdorff. Each member of We may now be viewed as a continuous function 
on X, since it has a well-defined value (zero) at [0] = [+oo] G X. 

Clearly, We_®R := {^p + c \ (p G Wg, c G M} forms a subalgebra of C{X) 
that contains the constant functions, and now separates the points of X (since the 
offending pair in [0, +oo] have been merged). We may now apply Stone- Weierstrass, 
to conclude that We © M is dense in C (X) under the supremum-norm. 

Returning to [0,+oo], this means that any ip G C([0,+oo]) for which ip{0) = 
(p{+oo) may be uniformly approximated by functions from H^£©]R, whence follows 
that = Co,o(M+). 

To finally prove the required property, let ip G Co,o(K+), and let e > 0. By the 
density of We, there is ^ G We such that \\(p — ip\\oo < e. If 9p = 0, we are done, 
since in this case. We = Yg. If 9p > 0, then ip may not satisfy the conditions at t = 
required for membership of the subset Ye, so we need to perturb Lp to achieve this. 
Place 

e := {w{0+) - ^(0), . . . ,^(^-+^^-2)(0+) - ^(^-+^--2)(0)) = (0, . . . ,0,e.„-i, . . . ,e.»+..-2) 

the latter obtaining since ^^■'-'(0) = for j < 0^ — 2 and also if Xe 7^ we must have 
^'^•'^(0+) = for such j. 

We now seek p(-) G We such that ip + p E Ye and \\p\\oo is suitably small. We 
slightly modify a procedure used in \7]. Try p(t) = f^poit) where m = 6^^ — 1, and 
Po(^) = Yl'i=o(^i^~^^^^'''^^ 'with q := 9p — 1, and a > so large that snp^yQt"^e~"'^ < 
e/(l + 1^1). Then, for k such that m<k<m + q = 9^ + 9p-2, have p^^\0) = 
[k\/{k — m)\]pl^~"^\o) (as well as p^^^(O) = for k < m). For ^ + p to satisfy the 
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required constraint at t = 0, we need 

^{k m.)^g^ _ j-^^ _ ^-jiy'^ij^^ A; = m,m + l,...,m + g 

or, equivalently, 

' 1 1 . . . 1 
1 2 ... q + l 
1 22 ... + 

12'?... {q + iy 

Let 7 denote a suitable norm of the inverse of the above square matrix. Place 
K{q) = (7(g) + l)(g + 1) — note that it depends only on q. Solving the above for the 
coefficients q yields po satisfying |po(^)| < K{q)\C^\e~'^^ for all t > 0, so that 

\p{t)\ < K(g)re~"*|e| < K{q)e for alH > . 

Thus it follows for this choice of p that (p + p &Yq, and 

ll^ + P- V^lloo < e+ IIpIIoo < {l + K{q))e. 

Since K does not depend on the choice of or ^ etc, it follows from the arbitrariness 
of e > that Ye, and hence Xg, is dense in Co,o(IR+), as claimed. 



Co 




OlCJml 


Cl 










q\U+J{{-l)%rn + q)\a'i) 



7.0.2 Case of = 1 

The same reasoning apphes, but without recourse to quotient spaces of [0, +00]. 

When p > 0, we obtain that Yg is dense in Co,a(M+), with a = w{0+). 

If p = 0, the constraints at t = in the definitions of Yg and Wg are absent, and 
arguing as above shows that Yg is dense in the entire space Co(IR+). 

This completes the proof of Proposition 13.21 



8 Appendix C: Proof of Lemma 14.3 

We start with a lemma. 

Lemma 8.1 Let /i G M(M+). Then, 

||/i + A5|| = ||/i|| + |/i({0}) + A|-|/i({0})|. 
In particular, z//i({0}) = (for example, when E U Q)V (BW ) we have \\fi + X6\\ 
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Proof: By definition, 



Wfi + \S\\ = l/i + A(5|(]R+) = sup < \fi + \S{I)\ : n partitions M+ into intervals > 

I /en J 

= li|!^o( ^ IM/)l + IM/o) + A| 

where /q denotes the unique member of 11 that contains 0. Since Jq i {0} for a 
subsequence of partitions, we have //(/q) — )■ //({O}), so that 



l/^ + A5|| = | 



hm (5^|M/)| + IM/o) + A|-|M/o)| ) 



= ||/i|| + |/i({0}) + A|-|/x({0})| 

□ 

Now, consider f = fos- Then, 

/#(/.) =rnax[aA + A5)] 

= max{aA | yU + A5 < 0, + A5|| < 1} 
= amax{A \ ^ + X6 < 0, \\^\\ + |A| < 1} by Lemma [8? 

Let n e U ®V ®W. If > 1, then + XS\\ = \\fi\\ + |A| > 1 for each A, so 
fosif^ + X6) = +00 for all A, and ffg^ii) = — oo = —f*g{fi). If /i is not a negative 
measure, then /i(-E') > for a Borel set E' in M (with ^ E), implying for all 
A G M, that (yU + X6){E) = fi{E) > 0, so fi + X6 also not negative, and again, 
/os(/" + A5) = +00 for such A, yielding = -cx) = 

If fl < and < 1 (so /i G dom/*^), then, since ^ + X6 < implies A < 
(take E = {0}), it follows that f^ifJ^) = 0, and hence equals —f*g{fi). 

This proves that /* = — /*^ on f/ © V" © W. The proof that /* = —fus is similar. 

To treat //;, we note an elementary lemma. 

Lemma 8.2 Let n E M(R+) with /i({0}) = 0. Then, (/x + X5)± = n± + X±5 for 
any A G M. 

Then, by similar arguments to those above, f^{fi) = —oc if /i ^ dom/^ = {/z | 
/i+(R+) < 1/2 & /i_(M+) > -1/2}, and for G dom/*, /#(/i) = a(l/2 - /i+(M+)). 
Thus the required relation for /a is also established. 

Similar reasoning applies for the remaining functionals. 
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